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Abstract 

We provide an analytical estimate of the effect of a spherical inhomogeneity on light 
beams that travel through it. We model the interior of the inhomogeneity through the 
Lemaitre-Tolman-Bondi metric. We assume that the beam source is located outside the 
inhomogeneity. We study the relative deviations of travelling time, redshift, beam area 
and luminosity distance from their values in a homogeneous cosmology. They depend on 
the ratio H = Hr^ of the radius ro of the inhomogeneity to the horizon distance 1/H. 
For an observer located at the center, the deviations are of order H 2 . For an observer 
outside the inhomogeneity, the deviations of crossing time and redshift are of order H 3 . 
The deviations of beam area and luminosity distance are of order H 2 . However, when 
averaged over all possible locations of the observer outside the inhomogeneity, they also 
become of order H 3 . We discuss the implications for the possibility of attributing the 
observed cosmological acceleration to the emergence of large-scale structure. 
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Introduction: The cause of the perceived acceleration of the present cosmological 
expansion has not been identified yet. An interesting possibility, that does not require 
the introduction of new ingredients to Standard Cosmology, is that the growth of inho- 
mogeneities in the matter distribution affects the astrophysical observations similarly to 
accelerated expansion in a homogeneous background. In particular, the luminosity dis- 
tance of faraway sources may be increased because of the propagation of light through 
inhomogeneous regions before reaching the observer. 

An unambiguous way to examine this possibility is through the study of the trans- 
mission of light in an exact inhomogeneous background. The analytical modelling of the 
Universe can only be approximate, and depends on the scale of the assumed inhomo- 
geneities. At length scales above O(10) hr 1 Mpc the density contrast is at most of 0(1). 
A popular choice for the background is based on the Lemaitre-Tolman-Bondi (LTB) met- 
ric [1]. The background has spherical symmetry, but can be inhomogeneous along the 
radial direction. The metric can be matched to the Friedmann-Robertson-Walker (FRW) 
metric at a certain radius r§. There are two possible choices for the location of the 
observer, which are consistent with the isotropy of the Cosmic Microwave Background, 
a) He/she could be located in the interior of the inhomogeneity, near its center [2]. b) 
He/she could be located in the homogeneous region, with the light travelling across several 
inhomogeneities during its propagation from source to observer [3, 4, 5]. 

In both cases, the size of the inhomogeneity tq determines its effect on quantities such 
as redshift and source luminosity distance. The relevant quantity is the dimensionless 
ratio H = tqH of tq to the horizon distance 1/H. Consistency with observations requires 
that H be of O(10~ 2 ), even though values larger by an order of magnitude have also 
been advocated for the explanation of the supernova data [2]. In the following we use 
perturbation theory in H in order to determine the dependence of the photon redshift 
and source luminosity distance on H, for both possible locations of the observer. 

Gravitational background: The LTB metric can be written in the form 

fa* = _ dt 2 + ^(MO dr 2 + R 2 { ^ r)dn 2^ (1) 

where dti 1 is the metric of a two-sphere, the prime denotes differentiation with respect to 
r, and f(r) is an arbitrary function. The function R(t, r) describes the location of a shell 
of matter marked by r at the time t. Through an appropriate rescaling it can be chosen 
to satisfy R(0, r) = r. 

The Einstein equations reduce to 

M'(r) = 47rR 2 p(t,r)R', (3) 

where the dot denotes differentiation with respect to t, and G = (l67rM 2 ) 1 . The gen- 
eralized mass function M(r) of the pressureless fluid with energy density p(t,r) can be 
chosen arbitrarily. 

We parametrize the energy density at some arbitrary initial time tj = as Pi(r) = 
p(0, r) = (1 + e(r)) po,i- The initial energy density of the homogeneous background is po,i- 
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If the size of the inhomogeneity is ro, the matching with the homogeneous metric in the 
exterior requires 4ir J r ° r 2 e{r)dr = 0, so that M.{tq) = 47rroPo,j/3. As we assume that 
the homogeneous metric is fiat, we also have /(ro) = 0. Discontinuities in f'(r) result in 
discontinuities in the derivatives of the metric functions. 

In our modelling we assume that at the initial time ti = the expansion rate Hi = 
R/R = R'/R' is given for all r by the standard expression in homogeneous cosmology: 
H 2 = p 0ti /(6M 2 ). Then, eq. (2) with R(0,r) = r implies that 

fir) = ^r 2 (l - 3M{r) ) (4) 

For our choice of f(r), overdense regions have positive spatial curvature and tend to con- 
tract, while underdense ones negative curvature and expand faster than the average. This 
is very similar to the initial condition considered in the model of spherical collapse. Even 
though we work with the particular choice (4) for f(r), we expect that our conclusions 
are valid for other variations of the LTB metric as well. These may include an arbitrary 
function io( r ) resulting from the integration of eq. (2). As this function appears in the 
combination t — to(r), it becomes irrelevant for large times. Also, the radial coordinate r 
is often redefined so that pi is constant. As this is only a gauge choice, we do not expect 
it to affect the physical behaviour. The eventual collapse or fast expansion of a certain 
region would be determined by its spatial curvature, as in our model. 

Optical equations: The optical equations [6] can be written as [4] 

1 d 2 ^A 1 



VI dX 2 AM 2 ' 



0°) -* 2 (5) 



do 2 dy/A (k 3 ) 2 R 2 ( 3M(r) . 

TTV + ^f-JT-' 7 = TTP). ( P TZryT I ' ( 6 ) 



d\ y/A d\ AM 2 V 47T.R 3 

where A is the cross section of a light beam, A an affme parameter along the null trajectory 
and k l = dx l /dX. The shear a is important when the beam passes near regions in which 
the density exceeds the average one by several orders of magnitude. Within our modelling 
of large-scale structure, applicable for scales above O(10) h~ l Mpc, the average density 
contrast is not sufficiently large for the shear to become important [4]. 

We assume that, even for general backgrounds, the light emission near the source is 
not affected by the large-scale geometry. By choosing an affine parameter that is locally 
A = t in the vicinity of the source, we can set d\[A~jdX 



= yHs- The constant £l s can 

A— 

be identified with the solid angle spanned by a certain beam when the light is emitted by 



a point-like isotropic source. This expression, along with y A =0, provide the initial 

A=0 

conditions for the solution of eq. (5). 

In order to define the luminosity distance, we consider photons emitted within a solid 
angle Q s by an isotropic source with luminosity L. These photons are detected by an 
observer for whom the light beam has a cross-section A a . The redshift factor is 1 + z = 
uj s /oj = kg /k®, because the frequencies measured at the source and at the observation 
point are proportional to the values of k° at these points. The luminosity distance is 
Dl = (1 + z)y/ A /Q s , with A Q the beam area measured by the observer for a beam 
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emitted by the source within a solid angle £l s . The beam area can be calculated by 
solving eq. (5). 

It is convenient to switch to dimensionless variables. We define i = tHi, r = r/ro, 
R = R/ro, where Hf = po j i/(QM 2 ) is the initial homogeneous expansion rate and vq gives 
the size of the inhomogeneity in comoving coordinates. The evolution equation becomes 

r_ 3M(f) m m 

R 2 4vri? 3 R 2 ' { ' 

with M = M/(pQ : irl) and / = 6M 2 / / (po.^o) = f/ H f, H i = i^r . The dot now denotes 
a derivative with respect to t. We take the affine parameter A to have the dimension of 
time and we define the dimensionless variables A = HiX, k° = k°, k 1 = k 1 /Hi, &r = r$k 3 . 
The geodesic equations maintain their form, with the various quantities replaced by barred 
ones, and the combination 1 + / replaced by H~ 2 + /. The optical equations take the 
form 

1 M = -W) 2 --* 2 ^ 



V 7 ! d\ 2 
2 dyfi 
dX ' ^[A d\ 



da 2 dVA _ 3/ 73 \2_ 5 /_ 3M \ 

H 7= Tf CT = 



with p = p/po,i and a = a/Hi. The initial conditions become dVA/dX _ = \fsT s /Hi = 

^/^^ s and VA =0, with A = H 2 A and Cl = H 2 fl. 

A=o 1 1 

The effect of the inhomogeneity on the characteristics of the light beam can be cal- 
culated analytically for perturbations with size much smaller than the distance to the 
horizon. These have Hi <C 1. In the following we use Hi as a small parameter in a pertur- 
bative calculation of the luminosity distance and redshift. For small inhomogeneities, the 
variation of the Hubble parameter during the crossing by the light beam is very small. As 
a result Hi is almost identical with the value H at the time of detection of the beam. We 
consider beams with /c 3 = that pass through the center of the spherical inhomogeneity. 
Beams with k 3 / can also be considered along the same lines, even though the calcula- 
tion is much more involved. 



Travelling time and redshift: The travelling time for a beam that propagates 
across the inhomogeneity has been calculated in ref. [4] up to 0{H 2 ). We denote by f s 
the location of the source and by t s the emission time of the beam. The travelling time is 

t- i s = ±H (R(i s , f a ) - R(i s , f)) + H 2 p R'{t s ,f)R{t s ,f)df 

-Hf (R(i s ,f s ) - R(t s ,f)) R(t s ,f) + O(Hf) (10) 

for incoming and outgoing beams, respectively. The leading term in the above expression, 
of 0{Hi), is the standard Doppler shift. It is non-zero whenever the observer has a peculiar 
velocity relative to a source in the homogeneous region. 

We can make a comparison with the propagation of light in a FRW background. In 
this case we have R(t, f) = a(t)f = R(i, l)f. Let us consider light signals emitted at f s = 1 
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and observed at the center (f G = 0) of the inhomogeneity. The peculiar velocity of such 
an observer is zero and the term of O(Hi) vanishes. The difference in propagation time 
within the LTB and FRW backgrounds is 

to - (Qfrw = H? t R'(i s ,f)R(t s ,f)dr - ^-R(i s , l)R(i s , 1) + O(tif). (11) 

For signals originating at f s = and detected at r a = 1 the time difference has the 
opposite sign. As a result, the time difference for signals that cross the inhomogeneity is 
of O(Hf). 

A similar expression can be derived for the redshift of a light beam that passes through 
the center of the inhomogeneity. One finds [5] 

ln(l + z)= ±Hi (R(t 8 ,r 8 ) - R(i s , f)) 

+H? jT R{i a ,f) (R(i s , f a ) - R(i s ,f)) df + O(Hf) (12) 

for incoming and outgoing beams, respectively. 

For signals originating at f s = 1 and detected at f = the redshifts obey 

\n(-l±l—) = Hf f 1 R(t s ,f)(R(t s ,l)-R(t s ,f))df 

V 1 + ZFRW J JO 

-^-R(i s ,l)R(i s ,l) + 0(Hf). (13) 

For signals originating at f s = and detected at f = 1 the r.h.s. of the above equation 
has the opposite sign. As a result, the redshift difference for signals that cross the inho- 
mogeneity is of O(Hf ). 

Beam area: The beam area obeys the second-order differential equation (8), whose 
solution depends crucially on the initial conditions. In certain situations, the symmetry 
of the problem permits an exact solution. For example, for signals emitted from some 
point f s at a time i = i s and observed at f a = we have [7] 

VA = (l + z)R(i s ,r s )VCi. (14) 

Similarly, for signals emitted from the center f s = and observed at f at a time i we 
have _ _ 

V / I= R(t ,f )VTl. (15) 

However, for a signal that crosses the inhomogeneity we need to integrate eq. (8) from 
f = to f with initial conditions determined by the propagation from f s to f = 0. These 
include not only VA, but d\f~A/dr as well. An exact analytical solution is not possible in 
this case, and we have to resort to perturbation theory in Hi. We have checked that the 
expressions (14) and (15) are reproduced correctly by our results, up to second order in 
Hi. 
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The optical equations (8), (9) can be written in the form 
d 2 VI 1 dk 1 dVA 



= ~^ Itt1 (tt) (16) 




d fa\ 1 dk 1 a 2 dVA a 3 f Rk 3 \ (_ 3M \ 

^W) + ~^~d\¥ + ^/J~dF¥ ~ zyw) ( ' 

The first term in the r.h.s. of eq. (16) is of O (H 2 ) because p = 0(1) and k° /k 1 = 
dt/df = Hi dt I dr = O {Hi). The term in the r.h.s. of eq. (17) is also of O (Hf) because 
Rk 3 /k 1 = HiRd(j)/dr = O (Hi). As a result, the second term in the r.h.s. of eq. (16) 
is O (Hf) and, therefore, negligible. The shear plays no role, except for cases in which 
the light passes very close to an extremely dense concentration of mass. At the length 
scales that we are considering the energy density is smoothly distributed, and the shear 
can be neglected. As the first term in the r.h.s. of eq. (16) generates the deviations of the 
luminosity distance from its value in a homogeneous background, we expect the overall 
effect to be of O (Hf) . In the following we confirm this expectation through an explicit 
calculation, assuming a simplified form of the energy density. 

We consider beam trajectories that start at the boundary of the inhomogeneity, pass 
through its center and exit from the other side. These have k 3 = 0. We express dk\/d\ 
in eq. (16) using the geodesic equation [4], and omit the shear. As the FRW metric is 
special case of the LTB one, no change of coordinates is necessary. In this way we obtain 

d 2 VI ( itiM R" Hff \ dVA _ 3 R' 2 r 
dr 2 + { ± ^THjf + dr ~ 2 P 1 + H 2 f ' ( 8) 

where the positive sign in the second term corresponds to ingoing and the negative sign 
to outgoing geodesies. 
We use the expansion 

VI = \fW + Hi^W) + H 2 ^W) + O(Hf), (19) 

and calculate VAF> in each order of perturbation theory. The travelling time is given by 
eq. (10). We can set t s = so the geodesic inside the inhomogeneity is i = —Hi(f — 1) 
for ingoing, and t = Hi(f + 1) for outgoing geodesies. We treat t as an 0(Hi) quantity. 

Central underdensity: We identify the initial time in the background evolution 

with the time of light emission: ij = t s = 0. This implies that R'(0,r) = R'(0,f) = 1. 
•-/ _ 

Also R (0, r) = 1. The initial configuration that we consider has Pi(0, r) = for r < r± 
and pi(0, f) = 1/(1 — r 3 ) for f > f\. From (7) we can calculate various derivatives of R at 
i=0: 

R(t, f) = R(0, f) + iR(0, f) + 0(H 2 ) = 1 + tR(0, f) + O(Hf), (20) 
^(i,f) = ^R(0,f) + O(H 3 ). (21) 



For r > n we have 



R (0, f) = —. \ 1 . , R (0, f) = — , /V 1 . . (22) 

v ' 2r 3 (ff - 1) ' V ' r 4 (rf - l) V ; 
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For r < r\ both R (0, r) and R (0, r) are zero. For the initial configuration that we 
assume, R is a continuous function of f. However, R is discontinuous at f = f\ and 
r = 1, while i? has 5-function singularities at the same points. 

The initial conditions for the solution of eq. (18) for an ingoing beam can be taken 
v^4(l) = 0, dVA(l)/df = —1, without loss of generality. To zeroth order in Hi, eq. (18) 
becomes d 2 V A^ /df 2 = 0, with solution V (f ) = — (r — 1) for ingoing and V A(°) (f ) = 
r + 1 for outgoing beams. To first order in Hi, eq. (18) gives dVAW /df = -2, with 
solution — 2r + 1 for ingoing and V A^> (r) = r 2 + 2r + 1 for outgoing 

beams. These results are the same as for the case of a homogeneous background. 

The effect of the inhomogeneity appears in second order in Hi. We obtain 

= -|p(0,f)Vl(0), (23) 

with the upper sign corresponding to ingoing and the lower one to outgoing geodesies. As 

we have already mentioned, for r < r\ we have Pi(0, r) = R (0, r) = R (0, r) = 0. 

The above equation can be solved analytically through simple integration, with the 
values at the end of each interval determining the initial conditions for the next one. The 

only non-trivial point is that the <5-function singularities of R at r = r\ and r = 1 induce 
discontinuities in the values of dV /df at these points. These must be taken into 
account in a consistent calculation. The discontinuities can be easily determined through 
the integration of eq. (23) in an infinitesimal interval around each of these points. The 
remaining calculation is straightforward. It must be emphasized that the discontinuous 
density profiles that we are considering can be viewed as limiting cases of continuous 
ones, when the transition regions become infinitesimally thin. The integration of eq. 
(23) around the corresponding values of r picks up the leading contributions arising from 
the transition regions. Including these contributions is necessary in order to reproduce 
correctly the exact expressions (14), (15). 

For a photon beam that starts from the boundary at f = 1, travels through the center 
of an underdensity at f = 0, and exits at the diametrical point with f = 1, we find 

VA&(f = 0) = 1-1 J l + l (24) 
4 f{ + f i + 1 

and 

yfmif = 1) = 5- -5—? . (25) 

ff + fi + 1 

Putting everything together, we find that, when the photon exits the inhomogeneity 
at f = 1, 

yfAir = l) = 2 + AH i+ U- _ 2 3 — ) H 2 + O(Hf). (26) 

The expressions for a homogeneous universe are obtained by setting f\ = 0. The beam 
area and the luminosity distance are increased by the presence of the inhomogeneity 
(fi ± 0). 



7 



We also mention that, if the beam is emitted at r = 0, it exits the inhomogeneity with 
\/W){f = 1) = 1/4 and \fW) (f = 1) = 3/4, in agreement with eq. (15). 

Central overdensity: The initial configuration that we consider has f>i(0,f) = 1/ff 
for f < f\ and pi(0, f) = for f > f\. 
For r > f i we have 

#(0,f) = l Jf(0,f) = -1 (27) 

while for f < f i we have 

R(0,f) = ~, R(0,f) = 0. (28) 

zr 1 

The expressions for v 7 !^ and \fW) are the same as in the case of a central underdensity, 
as they are not affected by the inhomogeneity. For V A^ we find 

yfW){f = 0) = 1- ~ (29) 

and 

v / I(2T(f=l) = 5-4- (30) 
r 1 

Putting everything together, we find that, when the photon exits the inhomogeneity 
at f = 1, 

v / I(f = l) = 2 + 4^+ (5-^jHf + 0(Hf). (31) 

The expressions for a homogeneous universe are obtained by setting f\ = 1. In this 
case, the beam area and the luminosity distance are reduced by the presence of the 
inhomogeneity (f\ ^ 1). The singularity for f\ — > is an artifact of the perturbative 
expansion. Clearly, the expansion in Hi breaks down when the coefficient of Hf diverges. 

The increase of the beam area by a central underdensity with a certain f\ can always 
be compensated by the decrease because of an overdensity with a different value f[. If 
one requires that f\ and f[ be equal, the solution is f\ = f[ = 2 -1 / 3 . In this case, the 
central underdensity and its surrounding overdense shell, as well as the compensating 
central overdensity and its surrounding underdense shell, all have equ al vol umes. 

If the beam is emitted at f = 0, it exits the inhomogeneity with v / A( 2 )(f = 1) = 1/4 
and (f = 1) = 3/4, exactly as in the case of a central underdensity. 

Flux conservation: We have seen that, when a light beam crosses a certain inho- 
mogeneity, the deviations of the travelling time i Q and redshift z from their values in a 
homogeneous background are of O(Hf), while the deviation of A is of O(Hf). As a result, 
the effect on the luminosity distance is of O(Hf). This conclusion holds for any beam 
going through the inhomogeneity, even if the crossing is not central. The analytical esti- 
mate has been verified through the numerical solution of the optical equations [4, 5]. In 
particular, a central crossing of a void-like inhomogeneity (with a central underdensity) 
results in the increase of the luminosity distance by an amount of O(Hf) [4, 5]. This 



8 



result is in agreement with the analysis of ref. [8] , in which a sequence of central crossings 
is assumed during the propagation of light from source to observer. On the other hand, 
if the inhomogeneity is crossed through the overdense region near its surface a decrease 
of the luminosity distance by an amount of O(Hf) takes place [4, 5]. 

The conclusion that the redshift is affected by an amount of O(Hf) has a very im- 
portant implication. If the redshift is not altered significantly by the propagation in the 
inhomogeneous background, the conservation of the total flux requires that the average 
luminosity distance be the same as in the homogeneous case. The energy flux may be 
redistributed in various directions but the total flux must be the same as in the homo- 
geneous case [9, 10]. The maximal deviation from exact flux conservation is determined 
by the effect of the inhomogeneity on the redshift, which is of O(Sf ). As a result, even 
though the effect on the luminosity distance for a single crossing is of O(Hf), the maximal 
average effect for beams originating in the same source and crossing the inhomogeneity at 
various angles is of 0(Hf). In the case of an underdensity, the increase of the luminosity 
distance for central beam crossings is compensated by a reduction for beams that travel 
mainly through the peripheral overdense shell. The opposite happens in the case of a 
central overdensity. 

The above conclusion has been verified numerically in ref. [5], both for central un- 
derdensities and overdensities. An equivalent conclusion is that the maximal statistical 
effect for light signals received from randomly distributed sources in the sky should be of 
O(Hf ), similarly to the effect on the redshift. The statistical analysis of ref. [5] confirms 
this expectation. 

Conclusions: The effect of spherical inhomogeneities on light emitted by a distance 
source depends on H = r$H. For an observer located at the center of a spherical inho- 
mogeneity, the deviations of travelling time, redshift, beam area and luminosity distance 
from their values in a homogeneous background are of 0(H 2 ). The luminosity distance 
is increased by the presence of a central underdensity, while it is reduced by a central 
overdensity. The increase in the luminosity distance if the observer is located near the 
center of a large void can by employed for the explanation of the supernova data [2]. An 
increase of 0(10%), as required by the data, would imply the existence of a void with 
size of O(10 3 ) h~ l Mpc. Numerical factors can reduce the required size, depending on the 
details of the particular cosmological model employed [2]. However, a typical void with 
size of O(10) h' 1 Mpc leads to a negligible increase of the luminosity distance. 

If the observer is located at a random position within the homogeneous region, the 
beam can cross several inhomogeneities before its detection. Each crossing produces 
an effect of 0(H 3 ) for the travelling time and the redshift. For the beam area and the 
luminosity distance the effect is of 0{H 2 ). However, flux conservation implies that positive 
and negative contributions to the beam area cancel during multiple crossings. The size of 
the maximal average effect of each crossing on the beam area and luminosity distance is 
set by the effect on the redshift, which is of O(H^) [9, 10, 5]. Photons with redshift ~ 1 
pass through ~ (l/H)/ro = H -1 inhomogeneities before arrival, assuming that these are 
tightly packed. As a result, the expectation is that the maximal final effect for a random 
position of the observer is of 0(H 2 ) for all quantities. This conclusion is supported by 
the numerical analysis [4, 5]. 

We mention at this point that, even for a random position of the observer, there is a 
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bias in the residual effect on the luminosity distance for a limited sample of sources. The 
bias is towards increased values if the Universe is dominated by void-like configurations. 
We did not discuss this point in this letter, as we assumed that the data sample is large. 
A detailed study can be found in ref. [5], to which we refer the reader for the details. 

We conclude that the presence of spherical inhomogeneities does not influence suffi- 
ciently the propagation of light in order to provide an explanation for the supernova data, 
unless their size becomes comparable to the horizon distance. It is possible, however, 
that relaxing the assumption of spherical symmetry for the inhomogeneities may increase 
the influence of the local geometry on the beam characteristics and provide an effect at 
a lower order in H. The crucial question is whether the influence of the inhomogeneities 
on the redshift can become larger than the effect of 0(H 3 ) predicted by our model and 
the Rees-Sciama estimate [11]. The modelling of the Universe as an ensemble of inho- 
mogeneities, glued together by a homogeneous region (the "Swiss-cheese" model), may 
be too constraining. Photons that cross an inhomogeneity enter an evolving newtonian 
potential from a homogeneous region, to which they subsequently return. Within this 
modelling, the residual effect cannot be much larger than of 0(H S ). The elimination of 
the intermediate homogeneous region may be necessary in order to produce a larger effect. 
This possibility poses formidable technical difficulties, but merits further investigation. 
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